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ABSTRACT 



We report the first detection of pure rotational transitions of TiO and TIO2 at (sub-)millimeter wavelengths towards the red supergiant 
VYCMa. A rotational temperature, r rot , of about 250 K was derived for TiCV Although T yot was not well constrained for TiO, it is 
likely somewhat higher than that of Ti02- The detection of the Ti oxides confirms that they are formed in the circumstellar envelopes 
of cool oxygen-rich stars and may be the "seeds" of inorganic-dust formation, but alternative explanations for our observation of TiO 
and Ti02 in the cooler regions of the envelope cannot be ruled out at this time. The observations suggest that a significant fraction of 
the oxides is not converted to dust, but instead remains in the gas phase throughout the outflow. 

Key words. Astrochemistry; Stars: winds, outflows; Circumstellar matter, Stars: individual: VY CMa; Line: identification 
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1 Copious amounts of dust, omnipresent in cool circumstellar en- 
vironments, are produced by stars at the final stages of their evo- 
lution, yet the formation of dust by depletion of molecules from 
the gas phase is poorly understood. In low-mass stars, the for- 
mation of dust is thought to be most efficient in the asymptotic 
giant branch (AGB), and in massive stars during the red super- 
giant phase. There are two major types of dust produced by cool 
stars. Carbon stars (e.g., C-type AGB stars and RCrB stars), 
which have more carbon than oxygen in their circumstellar en- 
vironments, produce mainly carbo n dust which nucleates from 
abundant aromatic hydrocarbons (Sedlmavr & Dominik 1995, 
and references therein). In oxygen-rich circumstellar envelopes 
(e.g., M-type AGB stars and red supergiants), essentially all car- 
Don is locked into CO and the dust is inorganic. How such in- 
organic grains start to form remains elusive. After ruling out 
species containing the more abundant elements (Si, Fe, and Mg) 
and taking into account infrared (IR) observations of circum- 
stellar dust shells in co ol oxygen-rich giants and supergiants, 
iGail & Sedlmavn (Il998l) concluded on the basis of a chemi- 
cal equilibrium model that the formation of inorganic dust at 
higher temperatures entails molecules containing less abundant 
elements. Gail & Sedlmayr showed that at temperatures ex- 
pected at the onset of dust formation, ~1200K, the refractory 
molecules TiO and T1O2 are the first species in the chain of reac- 
tions that produce titanium oxide clusters which are the "seeds" 
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* Based on observations carried out with the Submillimeter Array, 
IRAM Plateau de Bure Interferometer, and data acquired within the 
ESO program 67.B-0504. 



needed to initiate further dust formation at t emperatures below 
the condensation temperature of TiO (see also lJeong et al.l20"03l) . 

Prior to this work, T1O2 had not been identified in an as- 
trophysical object either at optical or radio wavelengths, but 
TiO has been obse rved at visible wav elengths in M-type stars 
for many decades (Me rrill et al.lll962l and references therein). 
The electronic bands of TiO are so prominent in late-type stars, 
that the Morgan-Keenan spectroscopic classi fication is based in 
part o n the relative strengths of these bands (Morgan & Keenan 
fl973h . However, the electronic bands in these studies are pho- 
tospheric absorption features, which do not provide any insight 
into the warm circumstellar chemistry. Occasionally the opti- 
cal bands of TiO are observed in emission from the warm cir- 
cumstellar g as, but such cases are extremely rare and d ifficult 
to interpret dBarnbaum et al.Hl 9961 iKamiriski et alJl2010h . Per- 
haps the most extreme example of the importance of TiO as an 
opacity source in oxygen-rich stars at optical wavelengths are 
high-amplitude mira variables, for which molecular absorption 
of TiO ma y be responsible for the variation of optical light by up 
to 8 mag dReid & Goldstonll2002l) . 

Radio astronomers had long sought TiO in the interstel- 
lar medium and circumstellar environments dChurchwell et al.l 
1980: iMillar et al.lll987t [Briinken et ai1l2008l) . because obser- 
vations of its rotational lines might yield the temperature and 
abundances more directly than at optical wavelengths. Although 
there have been several attempts to observe TiO at radio wave- 
lengths (most with single antennas), all were unsuccessful. In 
the past few years it was realized that images obtained with sen- 
sitive interferometers with sub-arcsecond resolution might have 
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sufficient sensitivity to observe the pure rotational spectrum of 
TiO and possibly Ti02- 

The rotational spect rum of TiO was m easured in the labo- 
ratory many years ago (Na miki et al.lll988l) . and the essential 
laboratory spectroscopic measurements needed to identify T1O2 
at optical and radio wavelengths were reported more recently. 
The rotational spectrum of TiQ2 w as first observed in the lab- 
oratory when Briinkenetal. (2008) measured 13 transitions to 
high accuracy in the centimeter band. Briinken et al. determined 
the rotational and leading centrifugal distortion constants that al- 
low predictions of the astronomically important transitions in the 
millimeter and submillimeter bands which subseque nt measure- 
ments confirmed are accurate to 0.5 km s _1 or better (iKania et al.l 
1201 lb . Unfortunately, the prospect for identifying TiCh at optical 
wavelengths appears difficult, because the most promising band 
at ~5300 A may be obscured by a more intense band of TiO. 

The red supergiant V Y CMa appeared to be the best source 
to search for TiO and Ti02 at (sub-)millimeter wavelengths, be- 
cause optical emission bands of refractory oxides including ScO, 
YO, VO, and TiO have long been known to be present in this 
extreme starfl and the (sub-)millimeter spectrum of VY CMa is 
unusually rich for an oxygen-rich source. In recent spectral line 
surveys at millimeter wavelengths at the Arizona Radi o Observa- 
tory (e.g.. lZiurvs et al.l l2007; Tene nbaum et al.1 12010). the num- 
ber of detected molecules in VYCMa increased to 21 includ- 
ing metal oxides (e.g., AlO and PO) and other exotic molecules 
(e.g., AlOH and PN). Some of the newly discovered molecules 
are refractory species important for dust formation in oxygen- 
rich CSEs. 

In 2010, we made a wide band line survey of VYCMa in 
the 345-GHzband with the Submillimeter Array (SMA). Present 
in this survey are rotational lines of both TiO and TiOa- The 
identifications of TiO and Ti02 were confirmed in 2012 by in- 
dependent measurements with the IRAM Plateau de Bure Inter- 
ferometer (PdBI) of lines within the 220.7-224.3 GHz range. In 
Sect.[2] we briefly describe the SMA and PdBI observations; in 
Sect.[3] the line identification and properties of the TiO and Ti02 
emission are presented; in Sect.|4] the rotational-temperature di- 
agrams for both molecules are analyzed; in Sect. [5] the (sub- 
)millimter TiO emission is compared with optical observations 
of this molecule; and finally, in the Sect.|6l the observed prop- 
erties of the emission are discussed within the broad context of 
inorganic dust formation. 

2. Observations 

Rotational lines of the two simplest Ti oxides (TiO and Ti02) 
were detected in a 345 GHz-band spectral line survey with the 
smaS The results of the survey together with a detailed de- 
scription of the observations and data reduction will be described 
elsewhere. Here we provide a brief description of the observa- 
tions. 

The spectral line survey of V Y CMa covers the entire fre- 
quency range between 279.1 and 355.1GHz, except for very 
narrow 32 MHz wide gaps at every 2 GHz interval. The observa- 
tions were done between 17 January and 2 February 2010. The 
typical time spent in each 4 GHz wide band was about 3.9 h, and 

1 Ap proximate luminosity (scaled to a distance of 1.2 kpc. lZhang et alj 
2012) is 3-lO 5 L , effective temperature r e ff =3200-3500 K, and ra- 
dius ^ = 1500-3000^ (see Appendix in Hump hreys et all 120071 : 
IWittkowski et alj2012h . 

2 The raw survey data are available at 
http://www.cfa.harvard.edu/cgi-bin/sma/smaarch.pl 



the rms noise for spectra extracted with an aperture of 1 " was 
about 70 mJy per resolution element of 0.7 km s _I . The noise is 
fairly uniform over the entire range of the survey except for the 
region between 321-327 GHz, owing to a strong telluric absorp- 
tion band of water. 

The SMA was used in its extended configuration to 
produce a synthesized beam with a full width at half- 
maximum (FWHM) of about 0'.'9. The phase center was at 
a(2000)=07 h 22 m 58.332 s ,(5(2000)=-25 o 46'03.17'.'2. All neces- 
sary calibrations (gain, flux, bandpass) were secured during the 
observations. The data were calibrated in the MIR-IDL pack- 
ageB After a careful inspection of the flux calibration, we esti- 
mate that its typical uncertainty is about 15%. The continuum 
map of VYCMa (0.6 Jy) was used to self-calibrate the data in 
phase and to correct for some of the phase variations. Most of 
the further data processing, including map pr oduction from ca l- 
ibrated visibilities, was performed in Miriad (ISault etal.1119951) . 

To confirm our identification of TiO and Ti02 with the SMA, 
we made an independent observation of a 3.6GHz-wide band 
near 222.45 GHz with the PdBI on 22 April and 13 October 
2012. The frequency range (220.65-224.25 GHz) was chosen 
so as to maximize the number of unblended and relatively strong 
lines of TiO and Ti02 which could be observed simultaneously 
in the 1.35 mm band. The PdBI was used in the D configuration 
with six antennas. Owing to the low declination of VY CMa, the 
maximum elevation is only 19°. As a result, the observations 
were affected by antenna shadowing and the synthesized beam 
is very asymmetric (4'.'9x2'.'2, with a position angle of the ma- 
jor axis of 158°). The observations of VYCMa were interleaved 
every 22 min by short integrations on nearby phase and ampli- 
tude calibrators (first 0648-165 and 0646-306, then switching 
to 0648-165 and 0823-223 when the elevation of 0646-306 was 
too low). Pointing and focus were checked about every 65 min. 
The observations were made with the dual-polarization single- 
sideband receiver (tuned in the lower side band) and the WideX 
correlator, which provides a resolution of 2 MHz (correspond- 
ing to 2.7kms~ 1 at the observing frequency). The data were 
calibrated using the GILDAS0 software package and the final 
processing included self-calibration in phase on the strong con- 
tinuum of VY CMa. The thermal noise in the combined data set 
is ~8 mJy/beam per 2.7 km s _I bin. 

3. Line identification and measurements 

In all, 220 spectral lines were observed in the SMA survey. Ini- 
tially we were able to assign all except about 30 on the basis of 
the close coincidence between the astronomical frequency and 
laboratory rest frequencies (v rest ) tabulated in the standard cat- 
alogs (Cologne Databas e for Molecular Spectroscopy, CDMS, 
Miilleret al. 2005, 2001; Jet Propulsion Laboratory Catalogue, 
iPickett et al.lll998l) . Of the remaining unassigned features, lines 
of TiO were relatively easy to recognize in the spectrum because 
they were narrower than most others, while Ti02 was one of the 
last species considered in our iterative identification procedure. 

The full spectrum obtained with PdBI is shown and de- 
scribed in AppendixlAl All of the lines of TiO and Ti02 pre- 
dicted to be present in the spectrum are detected in these very 
sensitive observations. 

Below, we describe the identification and present the spectra 
of TiO and Ti02- Data from the SMA survey are presented first 
and are followed by a description of the PdBI data, as this is the 

3 http://cfa-www.harvard.edu/~cgi/mircook.html 

4 http://www.iram.fr/IRAMFR/GILDAS 
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chronological order that influenced the analysis. A search for 
isotopologues in all the available data is described in Sect. 13.51 
Some constraints on the emission angular sizes and position are 
given in Sect. 13. 61 The evidence in support of the identifications 
is summarized in Sect l3.7l 

3.1. TiO in the SMA survey 

The TiO radical has a X 3 A electronic ground state that is split 
into three fine structure components with energies relative to the 
lowest 3 Ai component of 140 K ( 3 A 2 ) and 280 K ( 3 A 3 ). The ro- 
tational spectrum of TiO had been measu red to high accuracy 
in the laboratory by iNamiki et al.l (Il988|) and the dipole mo - 
ment was measured recently (^=3.34D: lsteimle~ & Virgo 2003). 
Starting with the 7=3— > 2 transition, roughly every 32 GHz there 
is a rotational transition (7'— »7) separated by 3-4- GHz from the 
corresponding transition in the next lower or higher fine struc- 
ture component. Eight transitions of TiO were covered in the 
survey: all three components in 7=9— >8 and 10— >9, but only 
two components were covered in 7=11— > 10 because the third 
(at 355.6 GHz) lies slightly outside the upper end of the observed 
range. 

Six lines of TiO were detected above the 3cr level in the SMA 
survey. Maps of the integrated emission, confirm that all six lines 
are from a source centered (within uncertainties) on the star (see 
Sect. 13. 61 ). Of the two missing transitions, the line at 291.0 GHz 
is located close to a gap and the baseline in the region near the 
line is uncertain. If the baseline in this spectrum is removed in 
a least squares fit with a polynomial, the edge of the spectrum 
near the gap is slightly lower and the fluxes in the frequency 
range of the expected feature of TiO are greater than 3 times the 
local rms noise. However, in the following discussion we treat 
the line at 291.0 GHz as though it had not been detected, and the 
measurements cited here are for the spectrum with a linear base- 
line subtracted (see below). The line at 323.3 GHz (7=10— >9 in 
3 A 3 ) is located in the frequency region affected by the absorp- 
tion band of terrestrial water and the noise is much greater than 
the expected intensity. As a result, none of the lines of the 3 A 3 
component was definitely detected. 

The most intense line of TiO at 348.2GHz may be blended 
with very weak emission from a line of 34 S02 at 348.1 17 GHz. 
Approximately 30 lines of SO2 were observed in the survey, and 
from these a rotational temperature of 85 + 10 K was derived. On 
the assumption that T rot is the same for 34 S02, we estimate that 
the contribution from the possibly interfering line of 34 S02 is 
<5% of the measured fluxes of TiO. We cannot rule out a priori 
that the lines of TiO (and those of Ti02, see Sect. l3.31 > might be 
blended with species whose carriers have not been identified in 
the survey, nor can we exclude entirely the possibility of blends 
with weak lines from species that have been identified. However 
the contribution to the measured fluxes from the latter should not 
be significant, because of the small aperture used in our analysis 
(l"xl"). For example, in a molecule like SO2, which has in- 
tense broad lines throughout the entire spectrum, the emission is 
greatest for larger apertures and transition s between the lowest 
rotational energy levels fcf. lFuetal.ll2012l) . 

Listed in Table[T]are the transitions of TiO which were cov- 
ered with SMA and PdBI, their peak and integrated fluxes, the 
rms of the noise in the region near the lines, and the detection 
levels (signal-to-noise ratios, S/N). The latter are expressed as 
the integrated flux density, /, divided by the uncertainty, A/, 
resulting from the lrms noise level, i.e., A/=rmsAV c h- yn^h, 
where AV c h is the channel width for which the rms is specified 
and n c h is the number of channels within the integrated inter- 



Table 1. Observed parameters of the TiO lines. 



Transition 1 * 


Vrest 


E u 


C b 

Jul 


Peak c 


Flux d 


rms e 


S/N* 


7'->7, Q 


(GHz) 


(K) 




(mJy) (Jy km s 


')(mJy) 




PdBI 


7->6, 1 


221.580 


41 


13.71 


74 


1.2 


5.4 


26.3 


7->6, 2 


224.139 


181 


12.85 


73 


1.3 


7.3 


20.5 


SMA 


9^8, 1 


284.875 


67 


17.77 


150 


1.8 


46 


8.5 


9-*8, 2 


288.156 


207 


17.11 


240 


2.3 


66 


7.6 


9^8, 3 


290.998 


354 


16.00 


<160 


<0.7 


50 


<3.0 


10-»9, 1 


316.519 


82 


19.80 


230 


2.4 


50 


10.6 


10-»9, 2 


320.159 


222 


19.20 


580 


3.1 


174 


3.9 


10->9, 3 


323.314 


369 


18.20 


<1470 


<14.0 


1021 


<3.0 


11->10, 1 


348.160 


99 


21.81 


440 


4.7 


73 


14.1 


11->10,2 


352.158 


239 


21.27 


290 


2.3 


125 


4.1 



Notes. (a) Q = 1, 2, 3 refers to the fine structure components 3 A 1; 3 A 2 , or 
3 A3. (b) The quantum mechanical line strength of the rotational tran- 
sition. <c) Peak flux density at a resolution of 2.7kms~' (PdBI) or 
0.8kms~' (SMA) for the full synthesized beam. (d) Flux integrated be- 
tween 6 and 32 km s -1 and within the full synthesized beam. (c) The rms 
per bin of 2.7 km s _I (PdBI) or 0.8 km s" 1 (SMA). <f) The signal-to-noise 
ratio. 



val. For those transitions that were not detected, the 3<x upper 
limit on the fluxes is indicated. Spectra of the individual lines 
are shown in Fig.Q] while Fig.|2]shows a combined profile of the 
lines that were detected with SMA, excluding the blended line at 
348. 16 GHz. The combined spectrum is an average of individual 
spectra weighted by the rms~ 2 and scaled by the reciprocal of the 
integrated flux. 

The line profile of TiO as seen in the SMA combined 
spectrum alone, extends over the range from about 3 to 
>55kms _1 (Pig.[2j. Most of the emission occurs between ~3 
and ~32kms~ 1 and is well reproduced by a Gaussian profile 
centered at V LSR =20+2 km s~' and with a FWHM= 1 8+2 km s" 1 . 
In addition, there is emission that forms a wing that is extended 
to the red of the main component. The redshifted emission is 
not discernible in the spectra of individual lines of TiO observed 
with SMA (Fig.CD. 

3.2. TiO in the PdBI data 

Two lines of TiO were covered by the PdBI observation are 
both well detected: the 7=7— >6 transition in the 3 Ai and 3 A2 
spin components at 221.6GHz 224.1GHz. The red wing of 
the 221.6 GHz line may be very weakly contaminated by Ti02 
emission (see below) and the blue wing of the 224. 1 GHz line 
may be affected by vibrationally excited SO2 (unlikely; see Ap- 
pendixlAl. The spectra of all TiO lines are shown in Fig. Q] and 
their parameters are listed in Tablefl] The line at 221.6 GHz ap- 
pears more triangular in shape than that at 224.1 GHz, possibly 
owing to the much different excitation energies of 41 and 181 K. 

The average line profile of the PdBI lines is compared to 
the combined SMA profile in Fig. [2] The PdBI profile is slightly 
broader than that observed with SMA. In the more sensitive PdBI 
observations, not only the existence of the red wing is confirmed 
but also a blue wing is apparent. They form a broad "pedestal" 
extending from -15kms _1 to 65kms~ 1 (with an uncertainty of 
about 5kms _1 in each of the limits). Because the ratio of the 
central, Gaussian-like component to the broad pedestal is smaller 
in the PdBI data, the main component appears broader. 
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Fig. 1. Spectra of all transitions of TiO covered in the SMA survey and 
the PdBI spectrum. The shaded area is the range over which the spec- 
tra were integrated. Overlayed on the measured spectra are theoretical 
Gaussian profiles (in blue) plotted on the assumption of LTE at 1000 K. 
The ordinate is flux density in Jy. Spectra correspond to the aperture of 
l"xl" (SMA data) or to the full PdBI beam of 4'.'9x2'.'2 (PdBI data). 
The SMA spectra were Hanning smoothed. 
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Fig. 2. Two top panels: Spectra of TiO and TiCh obtained by combin- 
ing several lines from the SMA survey (blue line) or PdBI observations 
(filled grey histogram). They are compared to a combined profile of 
NaCl emission (third panel) and a strong line of SO2 (16o,i6 — > 15 115 ) 
from the SMA survey (bottom panel). Spectra were smoothed and the 
peak flux density of each was normalized to 1.0. The dashed (red) ver- 
tical line marks the line center of TiO at Vlsr =20 km s , and the dotted 
(green) vertical lines mark 14, 35.5, and 48.0 km s~ ] identified as peaks 
of velocity components in the SMA combined profile of Ti02- The out- 
ermost dash-dotted lines (cyan) show the full velocity range of the Ti02 
emission. 



3.3. 770 2 in the SMA data 

Titanium dioxide, T1O2, is an asymmetric top with C2 V sym- 
metry and a substantial dipole moment along the b inertial axis 
(/Lfft=6.33 D: IWang et alj |2009). Owing to the two equivalent off- 
axis 16 bosons, half of the rotational levels are missing. Be- 
cause it is a fairly heavy molecule, many transitions lie within the 
range of the survey. We found 20 unblended lines and three pairs 
of lines that form tight blends of Ti02 emission. The blend at 
337.2 GHz is the most intense feature of T1O2 in the survey and 
one that provided the first indication that we may have detected 
this elusive molecule. In the blends at 331.2 and 337.2 GHz, 
the flux from each line in the blended pair is nearly the same, 
because both arise from levels at nearly the same energy above 
ground and the line strengths are comparable. 
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Table[2] lists the transitions detected, their flux, and the rms 
noise. It also contains the uncertainties of the rest frequencies, 
which typically are below 0. 1 MHz, but in a few cases are larger 
than 1 MHz. Shown in Fig. [3] are the features with the highest 
S/N, including the blended features at 331.2 and 337.2 GHz. 
Four unblended lines with the highest S/N (at 285.9, 347.8, 
350.4, and 350.7 GHz) were combined to yield an average pro- 
file by the same procedure as that for TiO (Fig.|2)- The line 
profile of Ti02 is broad, with most features extending over the 
velocity range from about to 62kms _1 . The strong blend at 
337.2 GHz suggests, however, that the range may extend as far 
as -15kms _1 . At least three velocity components at about 14, 
35.5, and 48.0km s _1 are apparent in the combined spectrum. 

3.4. Ti0 2 in the PdBI data 

The Cologne Database for Molecular Spectroscopy lists eight 
lines of Ti02 in the frequency range covered by the PdBI ob- 
servations. They are marked in Fig. IA.ll where the length of 
the markers indicates the relative line intensities expected un- 
der local thermodynamic equilibrium (LTE) conditions at 300 K, 
which is approximately the gas rotational temperature derived 
from the survey data (Sect.|4). Four are present in the observed 
spectrum. The one at 221.1GHz is blended with emission of 
34 SC»2, which we estimate is weaker than the 34 S02 line at 
221.7 GHz; therefore, the feature at 221.1 GHz is dominated by 
Ti02 emission. Two other lines, at 221.41 and 221.45 GHz (i.e., 
70kms~ 1 apart) overlap slightly; and another two, at 220.9 and 
223.5 GHz, are not detected, but are localized in spectral regions 
free of emission and sensitive upper limits on their emission can 
be derived. The remaining lines, at 221.55 and 222.96 GHz, are 
expected to have peak intensities much below the rms noise and 
coincide with emission lines of TiO and SO2. The four detected 
lines are included in Table|2] and the unblended ones are shown 
in Fig. [3] The integrated flux of the 221.1 GHz line includes the 
contribution of the 34 S02 emission and therefore is given as an 
upper limit in Table|2] 

When the PdBI line profiles are compared to each other 
(Fig.[3]), the 220.8GHz line is centrally peaked at 30kms _1 , 
while the other two lines are dominated by blue-shifted compo- 
nents peaking between 5 and 10km s~ l . The profiles of the two 
latter lines are very similar. They have larger-than-average un- 
certainty in their rest frequencies and a relative shift of the order 
of one spectral channel (~2.7kms _1 ) cannot be excluded. The 
emission in these two lines can be traced to at least -12kms , 
which when combined with the positive velocity components ob- 
served in some of the SMA lines gives the full velocity range of 
Ti02 emission between -12 and 62kms~ 1 - see the comparison 
of the combined profile^ in Fig. [2] The PdBI data confirm that 
the broad emission of Ti02 consists of several emission com- 
ponents that have different intensity ratios in individual tran- 
sitions, suggesting that the components are characterized by a 
range of excitation temperatures. There is some correspondence 
(although not exact) in the positions of the main components 
seen in both combined profiles, e.g., both are dominated by the 
central double (U-shaped) peak. The PdBI spectra and measure- 
ments of TiO and Ti02 provided in TablesQ] and [2] and Figs.Q] 
and[3]correspond to the full synthesized beam. 



5 In the average PdBI profile, some emission outside the -12 — > 
62kms~' range is seen, because the two averaged lines are not well 
separated in velocity. 



3.5. Isotopologues 

We did not observe lines of the rare isotopologues of TiO or 
Ti02- Because the S/N of the main isotopic species is <26, non- 
detection of the rare isotopologues is consistent with the stan- 
dard (solar and meteoritic) elemental composition in which the 
main isotope 48 Ti is at least 9 times more abundant than all other 
isotopes of Ti. 



3.6. SMA spatial constraints on the TiO and T1O2 emission 

Maps of integrated emission of TiO and Ti02 are shown in 
Figs. E] and [5] They are compared to maps of continuum emis- 
sion at 35 1+2 GHz, and that of NaCl 7=26^25 and H 2 u 2 = 
2 J k,..k, -^>2.7, —> 61.6, which are both expected to be compact 
sources (iMilam et al] 120071: iMenten et alJl2006h . The center of 
the continuum emission is marked by a white cross in the figures. 

We measured the position and sizes of the TiO and Ti02 
emission by performing fits of a two-dimensional (elliptical) 
Gaussian to the emitting regions. If the synthesized beam is ap- 
proximated by a Gaussian, then a deconvolved sizes of the emis- 
sion can be calculated from the fits. For this purpose, we used 
an automated procedure imfit in Miriad. We found that for some 
TiO transitions the main emission component (between 10 and 
30 km s _1 ) is a point source, i.e., its intrinsic size is much smaller 
than the beam. For a broader integration interval (6-32 km s _1 ), 
the average deconvolved size is 0'.'5 FWHM. This value has how- 
ever a large uncertainty owing to the very low S/N of the lines. 
For the Ti02 emission analyzed within 0-62 km s -1 , the average 
(weighted by the uncertainty of the Gaussian fit) deconvolved 
size is O'.'6+O'.'l FWHM, in agreement with FWHM=0'.'5±0'.'l 
derived separately for the 337.2 GHz feature alone, which has 
the highest S/N among all the lines analyzed. 

Although the position of the emission for each individual 
TiO and Ti02 feature usually agrees to within the (large) errors 
with the central position of the continuum source, in the case of 
the TiO line at 316.5 GHz, the measured center is displaced by 
— 0'.'2 in RA or 5 times the formal error of the fit; similarly, the 
strong Ti02 emission at 337. GHz is displaced by -0'.'2 in RA or 
4 times the error. The error- weighted mean values of the TiO and 
Ti02 positions for all the analyzed features have a similar offset 
with respect to the continuum source: (Aor, A6)=(-0'.'15±0'.'06, 
0'.'12±0'.'08) for TiO and (-0'.'20+0'/05, 0'.'08±0'.'03) for Ti0 2 , 
where the 3cr errors are indicated. The shift appears to be real 
and is distinguishable in Figs. [4] and [5] A similar offset is mea- 
sured for other relatively compact molecular emission, e.g., that 
of vibrationally excited H2O (see also NaCl; Fig.|5]l. The nature 
of the offset is discussed in Sect. 16. II 

The SMA measurements in TablesQ] and [2] and spectra in 
Figs.QJOare for data extracted with a square aperture of 1 "xl" 
centered on the continuum peak. With the above constraints 
on the emission size and position, this aperture encloses more 
than 95% of the emission thereby minimizing the contribution 
of noise in the spectra. We repeated the entire data analysis pro- 
cedure including the possible shift in the extraction and obtained 
essentially the same results as those here. The PdBI measure- 
ments and spectra, correspond to the large synthesized beam of 
4'.'9x2'.'2, and therefore do not provide additional constraints on 
the size of the emission. 
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Fig. 3. Spectra of unblended lines of Ti02 with the the highest S/N, and two blended lines (blue). The velocity scales of the blended lines refer 
to the average rest frequency of the two components. Parts of the spectra that contain strong emission from other species were blanked for clarity. 
The lines found in the PdBI observations are shown in the top three panels on the left. The dotted (green) vertical lines mark the position of the 
three main components of the TiC>2 emission identified in the combined profile (cf. Fig.|2j. The ordinate is flux density in Jy and corresponds to 
the aperture of l"xl" in the SMA spectra or to the full beam size (4'.'9x2'.'2) in the PdBI data. The SMA spectra were Hanning smoothed. 



3. 7. Summary of evidence in support of the identification of 
TiO and Ti0 2 

The evidence that we have observed pure rotational lines of TiO 
and Ti02 in VYCMa is overwhelming. Following the initial 
identification of both molecules in the SMA survey, two new 
lines of TiO and four of T1O2 were observed at the predicted in- 
tensities in a blind observation with PdBI. The agreement of the 
astronomical frequencies (Vlsr) with laboratory measured fre- 
quencies, which are accurate to a small fraction of the line width, 
is excellent. There are no missing lines for either molecule, nor 
are there any lines that are anomalously intense. As shown in 
Fig-Hi the intensities of the eight lines of TiO between 221.58 
and 352.16 GHz that were detected and the two upper limits (at 
291.00 and 323.31GHz), agree to within the measurement un- 
certainties with those derived from a simple LTE analysis (see 
Sect.|4|. For Ti02, approximately 20 unblended lines in the 
SMA survey with intensities within a factor of a few of the most 
intense lines at 350.4 and 350.7 GHz, are present as predicted 
(Tabled . Maps of the integrated emission of TiO in Fig.|4]con- 
firm that all six lines observed with SMA are from a compact 



source (approximately) centered on the star, with an angular size 
that is much smaller than our 0'.'9 synthesized beam, as expected 
for a species predicted to be localized in the inner outflow. As 
the sample maps in Fig.|5]show, the emission from T1O2 is simi- 
larly compact as that of TiO, and those of NaCl and vibrationally 
excited H2O which were observed with the SMA under the same 
conditions and configuration as that of TiO and Ti02- 

Circumstantial evidence consistent with the identification of 
these two titanium oxides is provided by the line profiles, and 
an estimate of the likelihood of a misassignment owing to inter- 
fering background lines. Despite the fairly low S/N, the profile 
shapes of all the lines assigned to TiO and Ti02 are consistent 
(Figs. Q] and [3). Those of TiO are narrow and centrally peaked, 
while lines of Ti02 are broad and typically asymmetric with mul- 
tiple velocity components. Our confidence in the identification is 
further enhanced by the very low density of lines in our 77 GHz 
wide survey of approximately 1 line every 380 MHz. The aver- 
age separation between lines in V Y CMa is equivalent to about 
25 linewidths of TiO, or 5 linewidths of Ti02 — i.e., there is 
very little of the overlapping of features that frequently occurs in 
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Table 2. Parameters of detected TiC>2 lines. 
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S/N d 


J K' a ,K' c — * Jk„,K c 


(MHz) 


(K) 




(Jy kms -1 ) 


(mJy) 




PdBI 


235,19^234,20 


220812.1 


(0.2) 


230 


12.857 


1.48 


11.0 


10.4 


44,0— >33,1 


221095.7 


(0.0) 


25 


3.476 


<1.77 


8.3 


16.5 


325,27 — >324,28 


221406.6 


(0.5) 


424 


19.365 


0.91 


7.5 


9.4 


284,24 —^283,25 


221458.4 


(0.3) 


323 


14.873 


1.45 


7.5 


14.9 


SMA 




282196.8 


(0.0) 


50 


2.489 


1.08 


46.1 


3.3 


2, 1 1 it 1 — ^ 2X ) a on 

i ' 1 1,Z1 ''"U/U 


284371.9 


(0.1) 


153 


19.115 


3.51 


49.5 


10.1 


20o IS — > 19^ 17 


284951.3 


(0.0) 


161 


9.737 


2.25 


47.2 


6.8 


257,19-^256,20 


285859.7 


(0.4) 
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14.106 


2.73 


42.1 


9.2 


22-7 1 ^ — >22a 1 a 


287253.6 


(0.2) 


239 


11.902 


2.41 


55.5 


6.2 


21 7 1 ^ — >2 1 A 1 A 

/ , 1 J D, 1 

4-^0 'I'J 1A 


290830.5 
290858.8 


(0.2) 
(7.8) 


222 
753 


11.157 
28.013 


} 2.39 e 


52.5 


4.7 




293001.6 


(0.1) 
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9.733 


2.69 


59.5 


6.4 




294443.0 
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59.5 


4.0 
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46.7 
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193,n — * 1 82, 16 
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64.1 


4.6 


22i,2i — >2l2,20 


310554.7 


(0.0) 
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16.308 


3.39 


82.4 


5.8 


23 1,23 — >22o,22 


310782.7 


(0.0) 
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21.124 


5.03 


84.3 


8.5 


27g,2o — >27 7 ,2i 

268,18 — >267,i9 


331211.8 
331240.1 


(0.7) 
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} 3.71 


48.9 


8.1 


398,32 — >397,33 


331381.8 


(2.5) 
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1.33 


48.3 


3.9 


123,9— >ll2,10 


331599.6 


(0.0) 


68 
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1.14 


52.6 


3.1 


243,21 — »234,20 


336469.4 
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3.77 


65.4 


8.2 
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337196.1 
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54.2 


37.2 


16s,8 ^167,9 
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1.80 


81.1 


3.1 


12 4 ,8-^ll3,9 


345580.9 


(0.0) 


75 


3.995 


1.43 


53.7 


3.8 


242,22—^233,21 


347788.1 


(0.1) 
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4.08 


69.0 


8.4 


355,31 — >354,32 


348661.4 


(1.3) 
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1.68 
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3.6 


26(),26->25i,25 


350399.0 


(0.1) 
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8.19 


71.4 


16.3 


252,24 — »24i,23 


350707.9 


(3.7) 


229 
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6.73 


70.2 
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SMA deblended* 


278,20^277,21 


331212.8 


(0.7) 
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14.862 


1.87 


48.9 


4.6 


268,18— ^267,19 


331240.1 


(0.6) 


327 


14.125 


1.84 


48.9 


4.5 


251,25—^240,24 


337196.1 


(0.1) 
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23.128 


9.29 


54.2 


22.8 


24i,23 — >232,22 


337206.1 


(0.0) 


212 


18.345 


7.67 


54.2 


18.8 



Notes. (a) Rest frequency (and its error). <b) Flux integrated between and 62 km s 1 within the full synthesized beam. (c) The rms per bin of 
2.7 km s (PdBI) or 0.8 km s (SMA). (d) The signal-to-noise ratio. <c) Lines are not included in Fig.|7] <f) Parameters for the deblended lines (see 
Sect.fOl). 



sources which have a higher density of lines (e.g., IRC+10°216 
and hot cores). Therefore, the probability that the emission fea- 
tures assigned to all 8 lines of TiO, or all 27 lines of TiC>2 are due 
to chance coincidences with lines of other species is exceedingly 
low in VY CMa: < 10" 10 for TiO and < 10~ 18 for TiO 2 

The velocity ranges in which the emission of TiO and 
Ti02 is observed are well within those of other firmly identi- 
fied molecules in VYCMa. For example, in Fig. [2] we com- 
pare the newly discovered emission to a profile of NaCl com- 
bined of many lines detected in the survey, and to the strongest 
emission line of SO2 (at 283.46 GHz), all extracted for the 



6 For TiC>2, the probability is even lower as more than 27 lines are 
seen. Above, we have presented only lines that are not contaminated by 
any other emission, while there are many examples of TiC>2 lines which 
interfere with emission lines of other species; the presence of TiC>2 in 
many such cases is evident when the blended profile is compared to a 
sequence of other uncontaminated profiles of the species. 



same aperture l"xl". Lines of NaCl are among the narrow- 
est features observed in VY CMa (Mila m et al.l2007l) . while S0 2 
shows broad, multicomponent emission typical of many sulphur- 
bearing species in this object. The weak emission of the two ti- 
tanium oxides is located almost exactly within the velocity range 
of the SO2 profile, while the main component of the TiO emis- 
sion covers the same range as the NaCl profile. This correspon- 
dence is additional strong evidence in support of our identifica- 
tion. 



4. Analysis of rotational-temperature diagrams 

We made an attempt to constrain the physical conditions of the 
emitting gas by ana lysing the Boltzma nn (rotational tempera- 
ture) diagrams (e.g.. ISnvder et al.ll2005l) for both molecules un- 
der the assumption of LTE and optically thin emission. The ro- 
tational temperature (Trot) and total column density (AO were 
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Fig. 4. All but the bottom right panel: Maps of emission integrated in the 6-36 km s~' range for all the TiO lines covered in the SMA survey. 
The columns are ordered (left to right) by 3 A,, 3 A 2 , 3 A 3 ; and the rows (top to bottom) by J = 9 — > 8, 10 — » 9, and 11 — » 10. The bottom right 
panel: map of the continuum around 351 GHz. The beam at half power is shown in the bottom left comer of each map. The contours are drawn at 
1.5cr and 3<x (black), 4.5<x (red), 6cr and 9<x (cyan), 120o~ (white), where a represents the noise levels measured in the integrated-intensity maps. 
The white cross (+) marks the peak of the continuum emission. 



obtained from results of least-squares fits to the measurements 
weighed by their uncertainties. The fluxes of the PdBI data were 
scaled to correspond to the aperture of l"xl". 

TiO: Shown in Fig.|6]is a rotational temperature diagram of 
TiO derived from measurements in TableQ] Owing to the large 
uncertainties in the fluxes of the eight lines because of noise 
and systematic effects, we were unable to constrain the rota- 
tional temperature when the line at 29 1 .0 GHz was omitted in the 
analysis - i.e., the line representing the best fit parameters has a 
negative slope which is nonphysical under LTE conditions and 
inconsistent with the upper limit on the 29 1 .0 GHz line. As men- 



tioned in Sect. 13. II the line at 291.0GHz is close to the detection 
limit. When the upper limit on this emission was taken into ac- 
count in the fitting procedure, we obtained r rot =1010±870K and 
N(TiO)=(6.9±0.9)-10 14 cnr 2 . This solution is shown in Fig.© 
and also in Fig.[TJ where Gaussian profiles are overlaid on ob- 
served lines to illustrate the relative intensity scale at the derived 
temperature. The quoted errors were derived from the formal 
uncertainties of the fit. While the rotational temperature remains 
poorly constrained, the uncertainty in column density is small. 

The failure to constrain r, ot from the diagram may imply 
that the assumption of LTE is not valid for the TiO gas. In- 
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Fig. 5. The same as Fig.[4]but for the strongest emission features of Ti0 2 and intensity integrated in the range 0-62 km s _I (adequately wider for 
blends). For comparison, emission of NaCl 7=26— * 25 and H 2 u 2 =2 JK a ,K c —^2,i —*&i,6 is shown. The contours are at 5, 1.5a (black), and 10, 15, 
20, 25, 35, 120cr (cyan). See Table|2]for quantum numbers corresponding to the Ti0 2 lines shown. 



deed, VYCMa displays optical TiO emission whose presence 
implies non- LTE effects, most likel y related to the strong stellar 
radiation (cf. Kam inski et al.ll20T3h . Large scatter of points on 
a rotational temperature diagram and inability to constrain r iot 
are problems we have encountered for several other molecules 
observed in the SMA survey, including NaCl and SiS for which 
many lines (>15) are observed, all at higher S/N. Although, op- 
tical and infrared pumping may be important in the rotational 
excitation of TiO, the uncertainties in the current measurements 
are too large to allow any conclusive statements. 

Ti02: A rotational temperature diagram of Ti02 derived 
from the measurements in Table|2]is shown in Fig. [7] For those 
lines that are blended and arise from levels of similar energy, the 



flux of the individual line was estimated from the intensity of the 
observed feature divided by the Einstein coefficient (A„/) of each 
line (see the lower part of Table|2] labelled "SMA deblended"). 
The deblended lines are plotted in the diagram but were not in- 
cluded in the fit. The best fit to all detected lines (SMA and 
PdBI) gives r rot =255±24K and A^(TiO 2 )=(7.5+0.9)-10 14 cm- 1 . 
This solution is also consistent with numerous 5<x upper limits, 
some of which can be seen in the diagram. 

Although the fit is satisfactory (the confidence level is above 
99%), it corresponds to the average gas parameters over the en- 
tire velocity range between and 62kms _1 and omits the most 
blue-shifted component in the PdBI spectra. The different ve- 
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Fig. 6. Rotational temperature diagram of TiO. The upper limit on the 
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locity components contributing to the Ti02 emission may have 
different temperatures. 

In the above analysis, the data from SMA and PdBI were 
combined despite the possibility that the molecular emission is 
time variable. We note, however, that temperatures calculated 
for PdBI and SMA data separately are consistent within uncer- 
tainties (T™ I =200±60K and 7^=305+35 K, respectively), 
justifying our fit to the combined data. 

5. TiO at radio and optical wavelengths 

To compare the rotational profiles of TiO observed at radio and 
optical wavelengths, we used an optical spectrum of VY CMa 
obtained i n 2001 with the Very Large Telescope (VLT) and re- 
ported in Kami riski et al.l d20 1 3h . Refering to the TiO line list of 
G. J. Phillipd we searched for single-component and well sep- 
arated rotational lines, and chose two examples in the y (0,0) 
band. A line at the rest wavelength /Ji a b=7076.93A is mainly 
formed by the Q3(J-34) transition and is affected by the pres- 
ence of components from higher vibrational bands (at 9kms~ 1 
with respect to /li a b); moreover, the wings of the line are blended 
with those of other nearby rotational lines. Another line, at 
/ l lab =7082.64A, is formed mainly by P 3 (26) and g 3 (39) (and 



available on the web site of R. L. Kurucz http://kurucz.harvard.edu 



Ri(l 18) of the (4,3) band) and is much better separated from the 
neighboring rotational components; its red wing, however, has 
some contribution from the weak ^ 3 (59) line (at 7kms _1 with 
respect to the main components). The structure of the lines is 
illustrated in Fig. [8] where the profiles are also compared to the 
combined SMA profile. 

The profile of the optical and radio emission is very similar. 
The optical lines are however slightly narrower than the com- 
bined SMA profile if results of Gaussian fits performed within 
3-32 km s 1 are compared; the average FWHM of the two opti- 
cal lines shown in Fig.|8]is 12±3kms _1 . Some other rotational 
lines measured are as narrow as 9 km s _1 FWHM, but due to un- 
known level of the photospheric fluxes these widths may be un- 
derestimated. For the same reason, the presence of the broad 
pedestal seen in the radio spectra cannot be investigated. The 
peaks of the optical lines are at Vlsr=22+1 kms~' and there is 
no systematic shift of line centers with /. Taking into account the 
uncertainty in the wavelength calibration of 3rms=0.7kms _1 , 
this peak velocity is formally consistent with the location of 
the center of the com b ined S MA profile (VLSR=20+2kms~ 1 ). 
Similarly. IWallersteinl (1 19861) derived V LS R=24±2kms~ 1 from 
positions of the TiO band heads, which to within the errors is 
consistent with the peak position of the radio emission. How- 
ever, when the profiles are overplotted, the optical lines ap- 
pear to be systematically shifted by about 2kms~ I with re- 
spect to the radio profile (note that a compensating correction 
was applied in Fig. [8). Some red-shift of the optical emission 
with respect to the radio emission is expected owing to scatter- 
ing of the optical radiation on grains within the dusty n ebul a 
(Ivan Blerkom & van Blerkomll978l:lHumphrevs et al.l2005l) . As 
noted in iKamiriski et al.l (120131) . this red-shift of optical lines 
must be less than about 3 km s _1 and does not affect considerably 
the profiles for the line of sight towards the sta r; both effects are 
signif icant only in the outer parts of the nebula (Humphrey s et al.l 
l200l . 

Rotational temperatures derived from our measurements 
of the pure rotational spectrum and previous optical studies 
of Ti O are consistent, althou gh the uncertainties are high in 
both. iPhillips & DavTsI (Il987l) derived a rotational temperature 
of 600 K with an uncertainty of about 150 K from the rota- 
tional structure of the y' system at 6200 A (B 3 n-X 3 A), how- 
ever there are two possible sources of systema tic error in the 
optica l measurements. First, the analysis of IPhillips & Davisl 
(1 19871) requires a revision with a new TiO line list containing 
improved oscillator strengths (e.g., ISchwenkel 17998). The new 
data can change the shape and intensity of the simulated spec- 
trum of the emission (and the underlying photospheric absorp- 
tion), Jhvisjdejr^g^evised estimates of the rotational tempera- 
ture (lAllard et al.l 1200 0). Second, because the optical emission 
spect rum of VY CMa and the bands of TiO are time variable 
(e.g., Wal lersteinll 1 9861) . the excitation temperature of TiO may 
be time variable too. For example, rotational temperatures de- 
rived from the orange (0,0) band of the A 2 n~X 2 S system of 
ScO differed by more than a factor of two on two differe nt dates 
(7^=380 K and 820 K in 1962 and 1966, respectively; [HirbJJ 
119741) . 

The above comparison provides further strong evidence in 
support our identification of TiO radio emission. The similarity 
of the optical profiles and the main (central) component of the 
(sub-)millimeter emission suggests that they arise in the same 
circumstellar gas. 
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Fig. 8. Rotational profiles of TiO emission in the y system near 7080 A 
(panels first and second from the top) are compared to the combined 
SMA profile of TiO (third panel from the top). Bottom panel: the Ti 1 
line at 5193 A is overlaid on the combined SMA profile. The vertical 
lines mark positions of rotational components of the (0,0) band, unless 
otherwise indicated in the labels in the Figure. All the optical lines were 
shifted arbitrarily by -2 kins' 1 to facilitate direct comparison with the 
SMA radio profile. 



6. Discussion 

6.1. Localization and abundances of the TiO and T/O2 gas 

The displacement of the center of the molecular emission with 
respect to continuum peak (Sect. l3.6l l brings into question the 
actual position of the star in our self-calibrated maps. The con- 
tinuum source is resolved in the SMA observations at a very 
high S/N (~200 for the peak flux in a single 2GHz-band) and 
is elliptical in shape with 0'.'5x0'.'3 FWHM. Therefore the ori- 
gin of the radio continuum is certainly not entirely photospheric. 
Indee d, models of the sp ectral energy distribution of VYCMa 
(e.g.. lHarwit et al.l l2001) suggest that at the observed frequen- 
cies the total flux is almost entirely reprocessed dust radiation, 
with the direct photospheric flux constituting only a few percent 
of the total emission. Our measured flux of 0.6 Jy at 320 GHz is 
consistent with those models. Although both dust and molecular 
gas are likely to have a very inhomogeneous distribution around 
VY CMa even at lowest spatial scales, we believe it is the center 
of molecular emission that can be identified as the position of the 
star itself. In support to this interpretation comes the observation 
that emission of transitions from highly excited states, like that 
of H2O 1)2=2 Jk u ,k c =52,3 - *6 1,6 with £„=2955K, peaks at this 
location; since it requires a nearby energetic excitation source, 
e.g., energetic radiation or high-density gas of high temperature, 
it is very likely located close to the star. Additionally, to pro- 
duce centrally peaked profiles (like those of TiO and NaCl) and 
emission center displaced with respect to the kinematical center, 



the gas would have to move tangentially with almost no radial 
component, which is very unlikely in this source. Moreover, 
Mulleret al. (2007) noticed a very similar spatial shift between 
the centroid of the continuum emission at 1.3 mm and the po- 
sitions of the masers of SiO, H 2 0, and OH in VYCMa. The 
positional shift between the continuum and molecular emission 
will be discussed in more detail in a separate paper. In the fol- 
lowing discussion we assume the observed molecular emission 
is spatially centered on the star. 

That the Ti02 e mission reaches the outfl ow terminal velocity 
(V<x>=40-45 km s' 1 , Humphreys et al. 2005) implies that at least 
part of the outflow is already accelerated to Vco- In case of TiO, 
the broad pedestal can be traced to high expansion velocities but 
due to a very low S/N it is uncertain whether it indeed reaches 
Vco (Fig.[2]i. The main Gaussian-like component of the TiO emis- 
sion, however, is very likely formed in the inner outflow because 
it is centrally peaked and much narrower than 2 V^. Similar nar- 
row molecular emission lines i n VY CMa have be en attributed 
to the accelerating w i nd (e .g., iMilam et alJ l2007 | ). Fo llowing 
iTenenbaum & Ziurvsl ( 120091) and lKaminski et al] d2013l) . to bet- 
ter localize the TiO emission region we use the semi-empirical 
wind-acceleration model in VY C Ma constrained from multi- 
epoch observations of H2O masers dRichards et alJll998h . This 
model assumes linear wind acceleration within 75-440 mas from 
the star and does not take into account the inhomogenities of the 
wind at these scales. Ignoring the weak pedestal, the main TiO 
component has a full width of about 40±5 km s' 1 , implying that 
the emission region extends out to r=28±10R* (R*=2000R o ) 
or 0'.'18+0'.'07 at 1.2kpc and the emission corresponding to the 
broad pedestal is located at even larger distances from the star. 

With the average deconvolved size of the of Ti02 emis- 
sion of FWHM=0'.'6, more than 99% of the flux is enclosed 
within a region of the radius of r=3cr/2=3 V2 In 2 FWHM=0'.'4 
or r=50(±10)R*. Using this value and the size of r=28R+ 
for TiO, we can correct the column densities derived in 
Sect. [4] for the corresponding beam filling factors. In this 
way, we obtain A^(TiO)=(6.7±0.8)-10 15 cm^ and jV(Ti0 2 )= 
1. 8+0.2)- 10 1 5 cm' 1 For a mass-loss rate of lO _4 M yr _1 
Dane hi et a and an average expansion velocity of 

lOkms' 1 for TiO and 20kms _1 for Ti02, the molecular abun- 
dances with respect to the column of hydrogen (assumed to 
be in molecular form only) are [/i]=log(A^xio/^VH)=-8.6 and 
[/2]=log(A^Tio 2 /-VH)= _ 9.5, but these estimates are very approxi- 
mate. 



6.2. Ti oxides and the formation of inorganic dust 

Our interferometric observations of TiO and Ti02 allow us to 
critically examine present ideas about the onset of nucleation 
and the formation of inorganic dust in oxygen-rich objects with 
high mass-loss rates. In the chemical equilibrium model of 
Gail &Sedlmayr( 1998), TiO exists in the photospheres of cool 
stars because the bond dissociation energy is high. As a result it 
is observed in regions with effective temperatures T e g <4000 K 
as well as in the cooler material just above the photosphere. Ti02 
is predicted to form in the reacti on of TiO with H 2O which is 
abundant in these environments dRover et alJl2010|) . Owing to 
the lower bond dissociation energy of Ti02, it is predicted to be 
more abundant than TiO only at lower temperatures. At temper- 
atures above about 1400 K, TiO is the dominant carrier of gas 
phase titanium; below this temperature T1O2 is the most abun- 
dant carrier of titanium; and at HOOK only about 10% of all 
titanium-bearing molecules are TiO. At slightly lower tempera- 
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tures, TiC>2 nucleates to produce larger oxides and their clusters 
thereby initiating the formation of dust "seeds". 

Dust for mation in VY CMa occurs at a radius^of 0'.'05-0'/12 , 
or 6-15 R* dMonnier et al.ll2000t ILTSidaner & Le Bertreiri 996). 
The prediction of iGail & Sedlmavrl (1 19981) that TiO should be 
located primarily inside the dust formation zone is inconsistent 
with our estimate that the size of the TiO emission region of 
28 R*. Unfortunately, our estimate of the rotational temperature 
(r rot =1010+870K, Sect.lXTT) is too uncertain to test the predic- 
tion of Gail & Sedlmayr that T ^;1 100K, but the relatively large 
observed line widths of TiO are inconsistent with the model. Al- 
though it is plausible that the outflow traced by the TiO gas is 
already accelerated to large velocities below the dust formation 
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it 



regioiQ it is more likely that the TiO emission does not origi- 
nate from inside the dust formation zone, but r ather from a much 
more e xtended region. Similarly, the model of IGail & Sedlmayr] 
(Il998l) is in poor agreement with the observed properties of 
Ti02- The temperature derived from the observations is about 
250 K, but is predicted to be about 4-5 times higher. The Ti02 
emission comes from a region (r^0'.'4), that is a few times larger 
than the size predicted for the dust formation region, although 
the uncertainties in both are quite large. Moreover, the offset of 
the peak emission for both molecules with respect to the center 
of continuum source cannot be easily reconciled with the predic- 
tions. 

Inconsistencies are also encountered when the observa- 
tions are referred to the equilibrium-chemistry model of 
ISharp & Huebneri fl990). In their model which includes con- 
densation, emission of TiO or Ti02 only arises from the warm 
gas (T >1400K) located in a thin layer above the photosphere. 
Simple equilibrium-chemistry models evidently cannot explain 
the observations of VY CMa. Indeed, molecular abundances cal- 
culated for molecular equilibrium have been demonstrated to fail 
observational tests i n AGB stars (ICherch neff 2006) an d for some 
speci es in VY CMa (Tenenbaum & Ziurvsl20 09; Kamihsk i et alJ 
120131) . A non-equilibrium approach is needed to explain the cir- 
cumstellar chemistry of supergiants but such models do not exist. 

Absorption bands of TiO are common features in optical 
spectra of cool photospheres. Although emission bands of 
TiO are rarely observed, they are present in VYCMa, and 
very rich TiO emission is seen in the spectra of UE qu and 
V4332 Sgr dBarnbaum et al.lll996LlKaminski et alJl2010h — two 
peculiar stars thought to have recently experienced a major cat- 
aclysm (a planet-engulfment and a stellar m erger, respectively; 
iGeballe et all 120051: iTvlenda & Sokerl 12 006]) . E mission of TiO 
was a lso identified in sever al B[e] supergiants (jZickgraf et al] 
1989) and a few protostars (Hillenbran d et al.ll2.012l) . including 
some that are known to be explosive objects. Recently, an emis- 
sion feature of TiO was discovered in the IR (around lOjum) 
in fou r S-type Mira variabl es which are large-amplitude pul- 
sators (Smolders et aDl2012h . A common property of these ob- 
jects (including VYCMa) is the relatively violent environment 
in which the TiO emission is observed. In VY CMa, it has been 
suggested that the outflow was formed by episod ic events of 
enhanced mass loss or "localized mass ejections" ( Smit h et alJ 



8 The dust formation radius is poorly defined owing to the asymmetry 
and probable variability of the star. Its location derived from observa- 
tions depends on the stellar parameters and dust type that were assumed, 
explaining in part the broad range of values cited in the literature. 

9 Outflows in luminous cool stars may be driv en by radiation pres- 
sure a bsorbed in molecular bands as proposed by Jorg ensen & Johnsonl 
d 19921 . and references therein), or the mass loss mig ht be driven by the 
convective/magnetic activity of the giant star (e.g., ISmith et ail 120011) 
which may lead to considerable turbulent motions very close to the star. 



2001 : Humphreys et al.l l2007) most likely involving shocks (cf. 
Kami hski et alJ 2013 ). If newly formed dust is destroyed close 
to the star by shocks, it might give rise to molecular emission 
at temperatures below ~1 100K. The observed variability of the 
optical TiO bands also suggests that the emission is linked to the 
stellar activity. 

However the kinematic separation of TiO and Ti02 implied 
by the line profiles in Fig.|2](i.e., TiO emission is highest at ve- 
locities where the Ti02 emission has local minima) would ap- 
pear to contradict this interpretation - one might expect that 
simple destruction of the dust would release TiO and Ti02 at 
the same location. A plausible explanation which is consistent 
with the observations is that when violent shocks release TiO, 
the Ti02 is destroyed but it reforms when water reacts with TiO 
farther from the star. 

Another possible explanation for the Ti oxides we observe is 
incomplete formation of the dust. The extended reflection neb- 
ula confirms that dust has been ef fectively formed in VY CMa 
for the last 1000 yr or more (e.g., iHumphrevs et al.| [2007). but 
owing to the high mass-loss rate (~10~ 4 M yr -1 ) a sufficient 
fraction of TiO and Ti02 might survive the nucleation pro- 
cess to account for the observed abundances in the cooler re- 
gion. The apparent stratification of TiO and Ti02 that we ob- 
se rve might then be unde rstood within the context of the model 
of IGail & Sedlmavrl d!998l). In addition, equilibrium-chemistry 
models of ISharp & Huebneri (119901) show that the titanium ox- 
ides remain relatively abundant throughout the outflow if no con- 
densation takes place. VYCMa may be an intermediate case 
between full versus no condensation. 

All of the models referred to here predict that titanium is con- 
verted to the oxides in the inner outflow, but emission features 
of atomic Ti are observed in the visible at the sys temic velocity 
(IWallersteinlll986t IWallerstein & Gonz alez 200lT). By examin- 
ing a pure emission feature of Ti 1 at 5193 A in the VLT spec- 
trum, we find that it covers the range from to 45kms~ 1 , i.e., 
is slightly broader than the central component of the radio TiO 
emission (see Fig. [8}. This implies that gas-phase Ti 1 extends 
farther from the star than the gas with highest TiO abundance. 
It has long been known that atom ic Ti is high ly depleted in the 
diffuse interstellar medium (e.g.. lStokeslll978l) . implying that Ti 
either forms molecules and remains in the gas phase, or it is ef- 
fectively locked into dust. The observed optical emission of Ti is 
further evidence that the condensation of Ti species in VY CMa 
is not as efficient as previously estimated, or that the Ti species 
extend to much larger radii than commonly assumed. Any future 
attempts to explain this anomaly should take into account those 
optical observations of atomic titanium. At this moment the ob- 
servational characteristics of Ti-bearing species in VYCMa re- 
main unexplained. 



6.3. Prospects 

The detection of TiO and Ti02 in the cooler region of the CSE of 
an oxygen-rich late-type star, implies that other small transition- 
metal bearing molecules such as VO, ScO, YO, and CrO might 
be found with sensitive interferometers in the (sub-)millimeter 
wave band. Observations of TiO and Ti02 with the Atacama 
Large Millimeter Array (ALMA) when the full sensitivity and 
angular resolution is attained, should resolve the region where 
they are present and establish whether there is evidence for strati- 
fication between species. This in turn would allow more detailed 
studies of titanium depletion to oxides and solids. 
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Kaminski et al.: Rotational spectra of TiO and T1O2 in VY CMa 
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Table A.l. Lines identified in the PdBI spectrum. 
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-8.9 


NaCl 
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Ti0 2 
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93 
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60 
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(S0 2 ) 
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23.0 
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/=29->28 
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22.4 


S0 2 
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9.9 
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64,2— > 73,5 
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35.0 


TiO 
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224.139 
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12.854 
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S0 2 


20 2 ,i8->19 3 ,i7 


224.265 
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Notes. (a) Quantum numbers characterizing the transition in the stan- 
dard notation. <b) The quantum mechanical line strength of the rotational 
transition. (c) Centroid position of the emission with respect to the rest 
frequency. 

Appendix A: The PdBI spectrum of VYCMa 

A spectrum extrac ted from the central pixel of the PdBI data 
is displayed in Fig JATFl The line identification for this spec- 
trum was based mainly on the results of the SMA su rvey, but 
al so on some archival SM A data (e.g., iFu et al.ll2012l) and that 
of iTenenbaum et all (1201 Ol) were consulted. The identified lines 
are shown in Fig. IA. H and listed in Table lA. ll In addition to the 
targeted TiO and Ti02 lines, the spectrum covers emission fea- 
tures of S0 2 , 34 S0 2 , NaCl, and KC1. Among these, KC1 had not 
been reported in VY CMa before (a full analysis of its emission 
will be presented elsewhere). Near 222.45 GHz, emission and 
absorption features that seem to compose a P-Cyg profile are 
seen. The emission peak can be tentatively identified with vi- 
brationally excited SO2 at V2-I, but the absorption occurs where 
two subunits of the WideX correlator are connected and is likely 
a spurious feature. 



10 The spectrum will be also available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 



http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-t-A/?? 
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Fig. A.l. The spectrum obtained with PdBI and corresponding to the central pixel of 0'.'65x0'.'65. The continuum was subtracted. The height of 
the lines marking TiC>2 positions (cyan), corresponds to the relative intensities for LTE at 300 K. The ordinate is flux density in Jy/beam. 
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